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Abstract — This document will review the fundamental ideas of linear algebra. 
We will learn about matrices, matrix operations, linear transformations and 
discuss both the theoretical and computational aspects of linear algebra. The 
tools of linear algebra open the gateway to the study of more advanced 
mathematics. A lot of knowledge buzz awaits you if you choose to follow the 
path of understanding, instead of trying to memorize a bunch of formulas. 

I. Introduction 
Linear algebra is the math of vectors and matrices. Let n be a positive 
integer and let R denote the set of real numbers, then R" is the set of all 
n-tuples of real numbers. A vector v G R" is an n-tuple of real numbers. 
The notation "G S" is read "element of 5*." For example, consider a vector 
that has three components: 

v = (vi,V2,V3) G (R,R,R)=R^ 

A matrix A G R'"^" is a rectangular array of real numbers with m rows 
and n columns. For example, a 3 x 2 matrix looks like this: 



ail 

£121 
£131 



£112 
£122 
£132 



The purpose of this document is to introduce you to the mathematical 
operations that we can perform on vectors and matrices and to give you a 
feel of the power of linear algebra. Many problems in science, business, 
and technology can be described in terms of vectors and matrices so it is 
important that you understand how to work with these. 

Prerequisites 

The only prerequisite for this tutorial is a basic understanding of high school 
math conceptaj like numbers, variables, equations, and the fundamental 
arithmetic operations on real numbers: addition (denoted +), subtraction 
(denoted — ), multiplication (denoted implicitly), and division (fractions). 

You should also be familiar with functions that take real numbers as 
inputs and give real numbers as outputs, / : R ^ R. Recall that, by 
definition, the inverse function f^^ undoes the effect of /. If you are 
given f{x) and you want to find x, you can use the inverse function as 
follows: [fix)) = X. For example, the function f{x) = h 



II. Definitions 

A. Vector operations 

We now define the math operations for vectors. The operations we can 
perform on vectors u = (ui,«2,U3) and v = {vi,V2,V3) are: addition, 
subtraction, scaling, norm (length), dot product, and cross product: 

U + V— {Ul + Vi,U2 + V2,Us + Vs) 
U~ V = {u-i- Vi,U2 ~ V2,U-i - Va) 

au — (aui, au2, auz) 



= \lul+ul+ul 



U - V — UiVi + U2V2 + UzV^ 
U X V —{U2VS — U3V2, U3V1 — U1U3, 111172 — U2V1) 

The dot product and the cross product of two vectors can also be described 
in terms of the angle 9 between the two vectors. The formula for the dot 
product of the vectors is u- v — ||m||||u|| cos^. We say two vectors u and 
V are orthogonal if the angle between them is 90°. The dot product of 
orthogonal vectors is zero: u - v — cos(90°) — 0. 

The norm of the cross product is given by ||m x i/H — \\u\\\\v\\ sin 9. The 
cross product is not commutative: u x v v x u, in fact ux v = —v x u. 



B. Matrix operations 

We denote by A the matrix as a whole and refer to its entries as aij. 
The mathematical operations defined for matrices are the following: 
• addition (denoted +) 



C = A + B 



subtraction (the inverse of addition) 
matrix product. The product of matrices A G 
is another matrix C G R 



and B G 



C = AB 



given by the formula 

n 

= aikbkj, 



ail 


£ll2 


0.21 


£l22 


£131 


£l32 



fell 

fe21 



612 
622 



fe=l 

£iiibii + ai2ti2i 
£i2ibii + a22b2i 

£l3lbll + £132^21 



aiibi2 + 012^22 

£121 bl2 + £I22&22 
(131 bl2 + £132622 



matrix inverse (denoted A ^) 
matrix transpose (denoted ^): 
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• matrix trace: Tr[^] = X^iLi 

• determinant (denoted det(y4) or \A\) 

Note that the matrix product is not a commutative operation: AB 7^ BA. 

C. Matrix-vector product 

The matrix-vector product is an important special case of the matrix- 
matrix product. The product of a 3 x 2 matrix A and the 2x1 column 
vector X results in a 3 x 1 vector y given by: 
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(ail, ai2) • X 

(a2i, 022) ■ X . (R) 
_(a3i, 032) • x_ 

There are twc[^ fundamentally different yet equivalent ways to interpret the 
matrix-vector product. In the column picture, (C), the multiplication of the 
matrix A by the vector x produces a linear combination of the columns 
of the matrix: y = Ax = xiAy.^i^ + X2A^.^2], where A[:,i] ^[■.,2] are 
the first and second columns of the matrix A. 

In the row picture, (R), multiplication of the matrix A by the vector x 
produces a column vector with coefficients equal to the dot products of 
rows of the matrix with the vector x. 

D. Linear transformations 

The matrix-vector product is used to define the notion of a linear 
transformation, which is one of the key notions in the study of linear 
algebra. Multiplication by a matrix A G R™^" can be thought of as 
computing a linear transformation Ta that takes n-vectors as inputs and 
produces m- vectors as outputs: 

: R" ^ R"". 



A good textbook to (re)learn high school math is minireference.com 



For more info see the video of Prof. Strang's MIT lecture: 'bit . ly/lOvmKcL 
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Instead of writing y = Ta{x) for the linear transformation Ta applied to 
the vector x, we simply write y — Ax. Applying the linear transformation 
Ta to the vector x corresponds to the product of the matrix A and the 
column vector x. We say Ta is represented by the matrix A. 

You can think of linear transformations as "vector functions" and describe 
their properties in analogy with the regular functions you are familiar with: 



function / : R ^ M 


<4> 


hnear transformation Ta : R" R™ 


input a; e R 




input s G R" 


output f[x) 


<^ 


output TA(f) = e R™- 


g° f^9{f{x)) 


<^ 


Tb{Ta{x)) = BAx 


function inverse /^^ 


<^ 


matrix inverse A~^ 


zeros of / 


<^ 


Af{A) = null space of A 


range of / 


<^ 


C{A) = column space of A = range of Ta 



Note that the combined effect of applying the transformation Ta followed 
by Tb on the input vector x is equivalent to the matrix product BAx. 

E. Fundamental vector spaces 

A vector space consists of a set of vectors and all linear combinations of 
these vectors. For example the vector space S = span {wi, {72} consists of 
all vectors of the form v — avi + I3v2, where a and P are real numbers. 
We now define three fundamental vector spaces associated with a matrix A. 

The column space of a matrix A is the set of vectors that can be produced 
as linear combinations of the columns of the matrix A: 

C{A) = {ye R'" \ y^ Ax for some f G R"} . 

The column space is the range of the linear transformation Ta (the set 
of possible outputs). You can convince yourself of this fact by reviewing 
the definition of the matrix-vector product in the column picture (C). The 
vector Ax contains x\ times the P' column of A, X2 times the 2°'' column 
of A, etc. Varying over all possible inputs x, we obtain all possible linear 
combinations of the columns of A, hence the name "column space." 

The null space N{A) of a matrix A G R'"^" consists of all the vectors 
that the matrix A sends to the zero vector: 

M{A) = {f G R" I Ax^Q}. 

The vectors in the null space are orthogonal to all the rows of the matrix. 
We can see this from the row picture (R): the output vectors is 0 if and 
only if the input vector x is orthogonal to all the rows of A. 

The row space of a matrix A, denoted TZ{A), is the set of linear 
combinations of the rows of A. The row space TZiA) is the orthogonal 
complement of the null space N(A). This means that for all vectors 
V G Ti-iA) and all vectors w G N{A), we have v ■ w = Q. Together, the 
null space and the row space form the domain of the transformation Ta, 
R" = N{A) ® TZ{A), where ® stands for orthogonal direct sum. 

F. Matrix inverse 

By definition, the inverse matrix A^^ undoes the effects of the matrix A. 
The cumulative effect of applying A^^ after A is the identity matrix 1: 

ri 01 

A-^A = 1 = 

0 



The identity matrix (ones on the diagonal and zeros everywhere else) 
corresponds to the identity transformation: T\{x) = lx = x, for all x. 

The matrix inverse is useful for solving matrix equations. Whenever we 
want to get rid of the matrix A in some matrix equation, we can "hit" A 
with its inverse A~^ to make it disappear. For example, to solve for the 
matrix X in the equation XA = B, multiply both sides of the equation 
by A-^ from the right: X = BA'^. To solve for X in ABCXD = E, 
multiply both sides of the equation by on the right and by A~^, B~^ 
and (in that order) from the left: X = C'^ B'^ A-^ ED-\ 



III. Computational linear algebra 

Okay, I hear what you are saying "Dude, enough with the theory talk, let's 
see some calculations." In this section we'll look at one of the fundamental 
algorithms of linear algebra called Gauss-Jordan elimination. 



A. Solving systems of equations 

Suppose we're asked to solve the following system of equations: 

Ixi + 2x2 = 5, 
3X1 + 9x2 = 21. 



(1) 



Without a knowledge of linear algebra, we could use substitution, elimina- 
tion, or subtraction to find the values of the two unknowns xi and X2. 

Gauss-Jordan elimination is a systematic procedure for solving systems 
of equations based the following row operations: 

a) Adding a multiple of one row to another row 

/3) Swapping two rows 

7) Multiplying a row by a constant 
These row operations allow us to simplify the system of equations without 
changing their solution. 

To illustrate the Gauss-Jordan elimination procedure, we'll now show the 
sequence of row operations required to solve the system of linear equations 
described above. We start by constructing an augmented matrix as follows: 
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The first column in the augmented matrix corresponds to the coefficients of 
the variable xi, the second column corresponds to the coefficients of X2, 
and the third column contains the constants from the right-hand side. 

The Gauss- Jordan elimination procedure consists of two phases. During 
the first phase, we proceed left-to-right by choosing a row with a leading 
one in the leftmost column (called a pivot) and systematically subtracting 
that row from all rows below it to get zeros below in the entire column. In 
the second phase, we start with the rightmost pivot and use it to eliminate 
all the numbers above it in the same column. Let's see this in action. 
1) The first step is to use the pivot in the first column to eliminate the 
variable xi in the second row. We do this by subtracting three times 
the first row from the second row, denoted R2 ^ R2 ~ SRi, 



\R2: 



3) We now start the backward phase and eliminate the second variable 
from the first row. We do this by subtracting two times the second 
row from the first row Ri ^ Ri — 2R2: 

1 0 1 
0 12 

The matrix is now in reduced row echelon form (RREF), which is its 
"simplest" form it could be in. The solutions are: xi = 1, X2 = 2. 

B. Systems of equations as matrix equations 

We will now discuss another approach for solving the system of 
equations. Using the definition of the matrix-vector product, we can express 
this system of equations l[TJ as a matrix equation: 



2) Next, we create a pivot in the second row using R2 
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This matrix equation had the form Ax = h, where A is a 2 x 2 matrix, x 
is the vector of unknowns, and 6 is a vector of constants. We can solve for 
X by multiplying both sides of the equation by the matrix inverse A~^: 



A ^Ax = lx = 


Xl 
X2 




A''b = 


[-1 ll 




■5' 
21 




T 
2 


But how did we know what the 


inverse matrix A ^ 


is 


? 







3 



IV. Computing the inverse of a matrix 
In this section we'll look at several different approaches for computing 
the inverse of a matrix. The matrix inverse is unique so no matter which 
method we use to find the inverse, we'll always obtain the same answer. 

A. Using row operations 

One approach for computing the inverse is to use the Gauss-Jordan 
elimination procedure. Start by creating an array containing the entries 
of the matrix A on the left side and the identity matrix on the right side: 
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2 


1 0 ■ 
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0 1 



Now we perform the Gauss-Jordan elimination procedure on this array. 
I) The first row operation is to subtract three times the first row from 
the second row: -R2 <— J?2 — 3-Ri. We obtain: 



IR2 
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-3 1 



2) The second row operation is divide the second row by 3: R2 
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3) The third row operation is Ri ^ Ri ~ 2R2 
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The array is now in reduced row echelon form (RREF). The inverse matrix 
appears on the right side of the array. 

Observe that the sequence of row operations we used to solve the specific 
system of equations in Ax = 6 in the previous section are the same as the 
row operations we used in this section to find the inverse matrix. Indeed, 
in both cases the combined effect of the three row operations is to "undo" 
the effects of A. The right side of the 2x4 array is simply a convenient 
way to record this sequence of operations and thus obtain A^^. 

B. Using elementary matrices 

Every row operation we perform on a matrix is equivalent to a left- 
multiplication by an elementary matrix. There are three types of elementary 
matrices in correspondence with the three types of row operations: 
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Let's revisit the row operations we used to find A~^ in the above section 
representing each row operation as an elementary matrix multiplication. 
I) The first row operation R2 ^ R2 — 3Ri corresponds to a multipli- 
cation by the elementary matrix Ei: 



C. Using a computer 

The last (and most practical) approach for finding the inverse of a matrix 
is to use a computer algebra system like the one at |live ■ sympy . org[ 

>>> A = Matrix ([[1,2], [3, 9]]) # define A 

[1, 2] 
[3, 9] 



>>> A.invO 

[ 3, -2/3] 
[-1, 1/3] 



# calls the inv metlnod on A 



You can use sympy to "check" your answers on homework problems. 

V. Other topics 
We'll now discuss a number of other important topics of linear algebra. 

A. Basis 

Intuitively, a basis is set of vectors that can be used as a coordinate system 
for a vector space. You are certainly familiar with the standard basis for 
the a::j/-plane that is made up of two orthogonal axes: the a;-axis and the 
j/-axis. A vector v can be described as a coordinate pair {vx,Vy) with 
respect to these axes, or equivalently as t; = v^i + Vyj, where i = (1, 0) 
and j = (0, 1) are unit vectors that point along the a::-axis and j/-axis 
respectively. However, other coordinate systems are also possible. 

Definition (Basis). A basis for a n-dimensional vector space S is any set 
of n linearly independent vectors that are part of S. 

Any set of two linearly independent vectors {61,62} can serve as a basis 
for ]R^. We can write any vector u G IR^ as a linear combination of these 
basis vectors v — V\e\ + ^262. 

Note the same vector v corresponds to different coordinate pairs depend- 
ing on the basis used: v = [v^, Vy) in the standard basis Bs = and 
V — {v\,V2) in the basis Be = {ei, 62}. Therefore, it is important to keep 
in mind the basis with respect to which the coefficients are taken, and if 
necessary specify the basis as a subscript, e.g., {v^,Vy)Bs "i" (^i)i'2)se- 

Converting a coordinate vector from the basis Be to the basis Bs is 
performed as a multiplication by a change of basis matrix: 



V 




1 






V 


<4> 


Vx 












Be 




Be 


Vy 





61 

61 





Vl 




V2_ 



■ 1 0' 




'1 2 




'1 2 




-3 1 




3 9 




0 3 


Mt 



|i?2 corresponds to a matrix E2 



EiA = 

2) The second row operation R2 

E2{EiA) 

3) The final step, Ri ^ Ri — 2R2, corresponds to the matrix E3: 

Ez{E2E^A) = 

Note that E-j,E2E\A = 1, so the product E-j,E2E-i must be equal io A' 
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Note the change of basis matrix is actually an identity transformation. The 
vector V remains unchanged — it is simply expressed with respect to a new 
coordinate system. The change of basis from the Bs-basis to the Be-basis 
is accomplished using the inverse matrix: — (s^[l]Be)^^- 

B. Matrix representations of linear transformations 

Bases play an important role in the representation of linear transforma- 
tions T: ]R" — > . To fully describe the matrix that corresponds to some 
linear transformation T, it is sufficient to know the effects of T to the n 
vectors of the standard basis for the input space. For a linear transformation 
T: 'M? — > ]R^, the matrix representation corresponds to 



T{i) T{j) 



As a first example, consider the transformation Tlx which projects vectors 
onto the s-axis. For any vector v — {vx,Vy), we have Ilx{v) = {vx,0)- 
The matrix representation of fl^: is 
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As a second example, let's find the matrix representation of Rg 
counterclockwise rotation by the angle 8: 
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cos H 
sin 6 
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cos 6 



The elementary matrix approach teaches us that every invertible matrix 
can be decomposed as the product of elementary matrices. Since we know 
A-^ = E3E2E1 then A = (A^^)"^ = (EaEiEi)-'^ = E'^E^^E^^^. 



The first column of Mug shows that Rg maps the vector i = IZO to the 
vector 1Z6 = (cos 6, sin^)^. The second column shows that Rg maps the 



vector J = 1Z|- to the vector lZ(v 



— (— sin 6', cos 9)^ . 
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C. Dimension and bases for vector spaces 

The dimension of a vector space is defined as the number of vectors 
in a basis for that vector space. Consider the following vector space 
5 = span{(l, 0, 0), (0, 1, 0), (1, 1, 0)}. Seeing that the space is described 
by three vectors, we might think that S is 3-dimensional. This is not the 
case, however, since the three vectors are not linearly independent so they 
don't form a basis for 5. Two vectors are sufficient to describe any vector 
in 5; we can write 5 = span{(l, 0, 0), (0, 1,0)}, and we see these two 
vectors are linearly independent so they form a basis and dim(5) = 2. 

There is a general procedure for finding a basis for a vector space. 
Suppose you are given a description of a vector space in terms of m vectors 
V = span{i;i, V2, ■ ■ ■ , Vm} and you are asked to find a basis for V and the 
dimension of V. To find a basis for V, you must find a set of linearly 
independent vectors that span V. We can use the Gauss-Jordan elimination 
procedure to accomplish this task. Write the vectors Vi as the rows of a 
matrix M. The vector space V corresponds to the row space of the matrix 
M. Next, use row operations to find the reduced row echelon form (RREF) 
of the matrix AI. Since row operations do not change the row space of the 
matrix, the row space of reduced row echelon form of the matrix M is the 
same as the row space of the original set of vectors. The nonzero rows in 
the RREF of the matrix form a basis for vector space V and the numbers 
of nonzero rows is the dimension of V. 



D. Row space, columns space, and rank of a matrix 

Recall the fundamental vector spaces for matrices that we defined in 
Section [lI-E[ the column space C{A), the null space M{A), and the row 
space TZ{A). A standard linear algebra exam question is to give you a 
certain matrix A and ask you to find the dimension and a basis for each 
of its fundamental spaces. 

In the previous section we described a procedure based on Gauss- Jordan 
elimination which can be used "distill" a set of linearly independent vectors 
which form a basis for the row space 7i{A). We will now illustrate this 
procedure with an example, and also show how to use the RREF of the 
matrix A to find bases for C[A) and N{A). 

Consider the following matrix and its reduced row echelon form: 



A 



13 3 3 

2 6 7 6 

3 9 9 10 



rref(A) 



13 0 0 
0 0 10 
0 0 0 1 



The reduced row echelon form of the matrix A contains three pivots. The 
locations of the pivots will play an important role in the following steps. 

The vectors {(1, 3, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1)} form a basis for n{A). 

To find a basis for the column space C{A) of the matrix A we need 
to find which of the columns of A are linearly independent. We can 
do this by identifying the columns which contain the leading ones in 
rref(yl). The corresponding columns in the original matrix form a basis 
for the column space of A. Looking at rref(A) we see the first, third, 
and fourth columns of the matrix are linearly independent so the vectors 
{(1,2, 3)"^, (3, 7, 9)"^, (3, 6, 10)"^} form a basis for C{A). 

Now let's find a basis for the null space, J^iA) = {a; G M'* j Ax — 0}. 
The second column does not contain a pivot, therefore it corresponds to a 
free variable, which we will denote s. We are looking for a vector with three 
unknowns and one free variable {x\, s,x-i,X4,Y that obeys the conditions: 
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Let's express the unknowns xi, xs, and X4 in terms of the free variable s. 
We immediately see that Xi = 0 and X4 = 0, and we can write xi = —3s. 
Therefore, any vector of the form (—3s, s, 0, 0), for any s G R, is in the 
null space of A. We write J^{A) = span{(-3, 1, 0, 0)"^}. 

Observe that the dim(C(A)) = dim(7?.(j4)) = 3, this is known as the 
rank of the matrix A. Also, dim(7^(^)) + dim{Af{A)) = 3 + 1 = 4, 
which is the dimension of the input space of the linear transformation Ta- 



E. Invertible matrix theorem 

There is an important distinction between matrices that are invertible and 
those that are not as formalized by the following theorem. 

Theorem. For an nxn matrix A, the following statements are equivalent: 

1) A is invertible 

2) The RREF of A is the nxn identity matrix 

3) The rank of the matrix is n 

4) The row space of A is R" 

5) The column space of A is R" 

6) A doesn't have a null space (only the zero vector N[A) = {0}j 

7) The determinant of A is nonzero AeX{A) 7^ 0 

For a given matrix A, the above statements are either all true or all false. 

An invertible matrix A corresponds to a linear transformation Ta which 
maps the rz-dimensional input vector space R" to the 71-dimensional output 
vector space R" such that there exists an inverse transformation T^^ that 
can faithfully undo the effects of Ta- 

On the other hand, an n x n matrix B that is not invertible maps the 
input vector space R" to a subspace C{B) C R" and has a nonempty null 
space. Once Tb sends a vector w G J^{B) to the zero vector, there is no 
Tg^ that can undo this operation. 

F. Determinants 

The determinant of a matrix, denoted det{A) or \A\, is a special way to 
combine the entries of a matrix that serves to check if a matrix is invertible 
or not. The determinant formulas for 2 x 2 and 3x3 matrices are 
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If the \A\=Q then A is not invertible. If |A| / 0 then A is invertible. 

G. Eigenvalues and eigenvectors 

The set of eigenvectors of a matrix is a special set of input vectors for 
which the action of the matrix is described as a simple scaling. When 
a matrix is multiplied by one of its eigenvectors the output is the same 
eigenvector multiplied by a constant Ae\ = \e\. The constant A is called 
an eigenvalue of A. 

To find the eigenvalue of a matrix we start from the eigenvalue equation 
Ae\ — Xe\, insert the identity 1, and rewrite it as a null-space problem: 

Aex = AleA (A- Al) ex = 0. 

This equation will have a solution whenever — Al| = 0. The eigenvalues 
of yl G R"*^", denoted {Ai, A2, . . . , A„}, are the roots of the characteristic 
polynomial p(A) — \A ~ Al|. The eigenvectors associated with the 
eigenvalue Ai are the vectors in the null space of the matrix (A — A,!). 

Certain matrices can be written entirely in terms of their eigenvectors 
and their eigenvalues. Consider the matrix A that has the eigenvalues of 
the matrix A on the diagonal, and the matrix Q constructed from the 
eigenvectors of A as columns: 
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then A = QAQ 



Matrices that can be written this way are called diagonalizable. 

The decomposition of a matrix into its eigenvalues and eigenvectors 
gives valuable insights into the properties of the matrix. Google's original 
PageRank algorithm for ranking webpages by "importance" can be 
formalized as an eigenvector calculation on the matrix of web hyperlinks. 

VI. Textbook plug 

If you're interested in learning more about linear algebra, you can check 
out my new book, the No BULLSHIT GUIDE TO LINEAR ALGEBRA. 
A pre-release version of the book is available here: gum. co/noBSLA| 



